arXiv:gr-qc/0503117v2 28 Jul 2005 


Linear stability of Einstein-Gauss-Bonnet static spacetimes 

Part I: tensor perturbations 

Gustavo Dotti and Reinaldo J. Gleiser 

Facultad de Matematica, Astronomia y Fisica, Universidad Nacional de Cordoba, 
Ciudad Universitaria, (5000) Cordoba, Argentina* 

We study the stability under linear perturbations of a class of static solutions of 
Einstein-Gauss-Bonnet gravity in D = n + 2 dimensions with spatial slices of the 
form E" x R + , E™ an n —manifold of constant curvature k. Linear perturbations 
for this class of space-times can be generally classified into tensor, vector and scalar 
types. The analysis in this paper is restricted to tensor perturbations. We show that 
the evolution equations for tensor perturbations can be cast in Schrodinger form, 
and obtain the exact potential. We use S —deformations to analyze the hamiltonian 
spectrum, and find an S-deformed potential that factors in a convenient way, allowing 
us to draw definite conclusions about stability in every case. It is found that there 
is a minimal mass for a D = 6 black hole with a positive curvature horizon to be 
stable. For any D, there is also a critical mass above which black holes with negative 
curvature horizons are unstable. 

PACS numbers: 04.50.+h,04.20.-q,04.70.-s 


I. INTRODUCTION 

The analysis of the properties and behavior of gravity in higher dimensions has become 
in recent years a major area of research, motivated in particular by developments in string 
theory. Among others, the Einstein-Gauss-Bonnet (EGB) gravity theory has been singled 
out as relevant to the low energy string limit [1], The EGB lagrangian is a linear combination 
of Euler densities continued from lower dimensions. It gives equations involving up to second 
order derivatives of the metric, and has the same degrees of freedom as ordinary Einstein 
theory. An appropriate choice of the coefficients in front of the Euler densities enlarges the 
local Lorentz symmetry to local ( A)dS symmetry [2, 3]. A number of solutions to the EGB 
equations, many of them relevant to the development of the AdS — CFT correspondence 
[4], are known, among them a variety of black holes in asymptotically Euclidean or ( A)dS 
spacetimes [5-9]. These were found mostly because they are highly symmetric. Analyzing 
their linear stability, however, confronts us with the complexity of the EGB equations, 
since the perturbative terms break the simplifying symmetries of the background metric. 
The linear stability under tensor perturbations of higher dimensional static black holes in 
Einstein gravity was studied in [10]; the stability of higher dimensional rotating Einstein 
black holes is analyzed in [11], The quasinormal modes of higher dimensional black holes 
are analyzed in [12] for Einstein gravity and in [13] for EGB gravity. 
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In this paper we consider spacetimes that admit locally a metric of the form 

ds 2 = — f(r)dt 2 + g(r)dr 2 + r 2 gijdx l dx\ ( 1 ) 

where gijdx l dx J " is the line element of an n— dimensional manifold £” of constant curvature 
ft = 1,0 or —1. Linear perturbations around (1) can be conveniently classified, following 
the scheme proposed in [14], into tensor, vector, and scalar perturbations. The ft = 1 case 
£" = S n gives, for appropriate / and g, cosmological solutions, as well as higher dimensional 
Schwarzchild black holes. The stability of these solutions under tensor perturbations was 
studied in [15]. In this paper we provide the details of the calculations leading to the 
results in [15] as we extend them to the cases ft = —1, 0. In Section II and Appendix A we 
introduce tensor perturbations around (1) and calculate the variation of the Riemann tensor, 
then in Section III we review the basics of Einstein-Gauss-Bonnet theory (EGB), exhibit 
known solutions of the form (1) from [5-8], obtain the perturbative equation for the tensor 
mode, and reduce it to a Schrodinger equation. In Section IV we classify the EGB solutions 
(1). A number of different possibilities arise depending on the spacetime dimension, the 
value of the cosmological constant and the strength of the coupling of the Gauss-Bonnet 
term (string-tension). Compact manifolds £ of negative (null) curvature can be obtained 
by taking quotients of hyperbolic space (Euclidean space) by appropriate discrete isometry 
groups, and black holes having such manifolds as horizons can be constructed in EGB gravity 
(for black holes with exotic horizons see, e.g., [8]). The stability of cosmologies and black 
hole solutions is studied in Section V using the S-deformation approach [16]. In spite of 
the complexity of the original Schrodinger potential, an S-deformed potential is found that 
factors in a convenient way and allows us to draw definite conclusions about stability in every 
case. Our preliminary work on vector and scalar perturbations [17] seems to indicate that 
this factorization is peculiar of the tensor mode. Conclusions about tensor perturbations 
can be found in Section VI. 


II. TENSOR PERTURBATIONS OF A CLASS OF STATIC SPACETIMES 


As stated in the previous Section, in this paper we consider spacetimes with metrics locally 
given by (1). We use a,b,c,d,... as generic indices, whereas are assumed to 

take values on A bar denotes tensors and operators on 
The non-zero Riemann tensor components of the metric (1) are: 
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The non-zero Ricci tensor components are 
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We study perturbations around (1) of the form, 

g a b —> g a b + h ab . (4) 

Indices of h ab are raised using the background metric, therefore 5g ah = —h ab . The first order 
variation of the Riemann tensors is: 


SR ab cd = - {Rab df h f c - R ab cf hf d + ( V b V c h a d - V a V c h b d ) + (V a VV - V b V d h a c ) } (5) 


For transverse ( y a h ab = 0) traceless ( g ab h ab = 0) perturbations (5) gives 

SRa C = \{-V d V d h a c - Rjh f c + RfhJ - 2R ad ^h f d } 
SR = -R d f h f d 


( 6 ) 

(7) 


from where 

SR a b = SR a c g C b + R a c h cb = — -^V d 'Vdh ab + - ( R a ^h b f + R b f h a f ) — R a kbfh' fk = L h) a b ( 8 ) 

A l being the Lichnerowicz operator. The transverse traceless condition does not restrict the 
perturbation, it (partially) fixes the gauge. Linear perturbations can be classified into tensor, 
vector, and scalar perturbations [14]. Tensor perturbations are specific of higher dimensional 
( D > 4) spacetimes and are the ones studied in this paper. They satisfy h ab = 0 unless 
(a, b) = ( i,j ). The non-zero components V a h bc for such a tensor are 

2 1 

^ t^ij d t hij V r - S,ij d r hij hij 71 % hj r h, j Sjk 71 j J ij ^, (9) 

Thus, tensor perturbations satisfy the conditions g^Vihjk = 0 and g l3 h t] = 0 (transverse 
traceless on £”) if and only if g ab W a h bc = 0 and g ab h ab = 0 (transverse traceless on the 
spacetime). Since transverse traceless tensors (TTT) on £” can be expanded using a basis 
of eigentensors of the Laplacian, we need only consider TTT perturbations of the form 

hij(t , r, x ) = r 2 0 (r, t)hij(x) ( 10 ) 

where r 2 is factored for later convenience, and 

V ; 'T/,//,, = 7 hij, V%j = 0, g t3 hij = 0 (11) 


Note that, since £” is a manifold of constant curvature, an eigentensor of the Laplacian on 
£” with eigenvalue 7 is also an eigentensor of A^, the Lichnerowicz operator on £” [10], 
with eigenvalue A given by 

A = 2 Kn — 7 (12) 

Solutions to equation (11) in the case £” = S n can be obtained from [18]. From equations 
(4)-(ll) we get the non trivial components of the variations of the Riemann tensor, the Ricci 
tensor and the Ricci scalar. These are displayed in Appendix A. 
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III. TENSOR PERTURBATIONS IN EGB GRAVITY 


The Einstein-Gauss-Bonnet (EGB) vacuum equations are 

0 = Q b a = AG(o) b a + G( i) fe a + aG( 2) b a (13) 

Here A is the cosmological constant, G( o) , = g ab the spacetime metric, Gp) , = R a b — \Rg ab 
the Einstein tensor and 


a 

G (2 ) b 


_ 75 de d 

— n c b h . 


'de 


2 R d c R, 


da 


cb 


2 R b c R c a + RR b 


-5 a b {R cd ef R ef cd - 4R d R d c + P 2 ) 


(14) 


the quadratic Gauss-Bonnet tensor. These are the first in a tower G( s ) b a , s — 0,1, 2, 3,... of 
tensors of order s in R ab cd given by Lovelock in [19]. As shown in [19], the most general 
rank two, divergence free symmetric tensor that can be constructed out of the metric and its 
first two derivatives in a spacetime of dimension d, is a linear combination of the G^ b a with 
2s < d. Here we consider the static spacetimes given by (1). These are foliated by spacclike 
hypersurfaces, orthogonal to the time-like Killing vector, that contain a submanifold of 
dimension n — D — 2 (D the spacetime dimension) of constant curvature k — 1,0 or —1. 
Inserting (2) in (13) we find that (1) solves the EGB equation (13) if [7] 

= fij) = k - r 2 ip(r) (15) 

and 'ijj(r) satisfies 
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From (3) and (15), the Ricci scalar for this solution is 

R = (n + 2)(n + l)V-(r) + 2r(n + 2 ) AU + 


dr 


dr 2 


(16) 


(17) 


TTT perturbations around this solution produce first order variations of the tensors 
G( s ) b a ,s = 0 , 1,2 which are trivial unless (a, b) = ( i,j ). Setting g — 1/f and using the 
equations in Appendix A gives 


and 
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For later simplicity, we introduce three functions Kj(r), defined by 


SG {2) / = | (0 - /V'J Ad + <j>'K 2 + 4>K 31 V. 
Perturbations around a solution of (13) satisfy the equation 

5G {1)a b + ahG (2)a 6 = 0 

which, after setting 4>(r,t ) = e u;t x(r) gives a second order ODE for 0(r) 

0 = -f 2 x"(r) +p(r)x'(r) + (q(r) + u 2 )x(r) 
2arfK 2 - rff - rtf 2 
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and switching to “tortoise” coordinate r*, defined by dr*/dr = 1//, this ODE can be cast 
in the Schrodinger form, 

d 2 <h 

- 2 + E(r(r*))$ = -cu 2 $ = A$ (28) 

dr* 

The spacetimes (1) will therefore be stable if (28) has no negative eigenvalues. On the other 
hand, properly normalized eigenfunction of (28) with suitable boundary conditions (see, e.g. 
[10] for details) having a negative eigenvalue (E < 0), signals the possibility of an instability. 
The explicit form of K(r ) is 


K (r) = r n / 2 1 yC 2 + a(n — 2 ) ^(n — 3)(« — /)— 


(29) 


The explicit form of the potential V (r) as a function of r and the parameters of the theory 
is rather lengthy. We notice however that the function q in (25) is, 


and the potential is given by, 


af")r 2 + a(n — 3) [(n — 4)(« — /) — 2r/' 
r 2 + a(n — 2 ) [(n — 3 )(k — f) — rf] 


V(r) 


, f d 

? + Jl dr 



(30) 


(31) 


V(r), given by (31) is the exact potential of the Schrodinger-like stability equation for the 
spacetime (1) in EGB gravity. This includes EGB blackholes of arbitrary mass and cos¬ 
mological constant, as well as cosmological solutions of the EGB equations that result by 
setting fi —— 0 in (16). It generalizes the k = 1 case first presented in [15], and it is readily 
seen to reproduce the potentials in [10] in the a = 0 (Einstein gravity) limit, a case that 
was extensively studied by Kodarna and Kodama and Ishibashi (see, e.g., [16] and references 
therein). The restricted cases in [20] and [21] can also be studied using (31). 
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IV. CLASSIFICATION OF MAXIMALLY SYMMETRIC STATIC SOLUTIONS 


A classification scheme for the solutions of the EGB equations is introduced below fol¬ 
lowing Whitt [7]. It should be kept in mind that a particular EGB theory is defined once 
the values of the space-time dimension n + 2 , the cosmological constant A and a (assumed 
different from zero), are given. A particular symmetric solution (1)-(15)-(16) of an EGB the¬ 
ory further requires the specification of the discrete index k and of the integration constant 
// in (16). Solutions are classified according to their singularities, horizons and asymptotic 
behaviors. To analyze singularities we rewrite the Ricci scalar (17) entirely in terms of 0. 
This is done using (16) and its first two derivatives together with (17). We arrive at 


R = { n{n — l)(n — 2) [n 2 (n + 3)(n + l)(n — l) 2 (n — 2) 2 a 2 0 4 

+4 n 2 (n + l)(n — l)(n — 2 )a0 3 + 8n(3 + 2 n)(n — 1 ){n — 2)Aa0 2 — 2 n 3 (n + 1)0 2 

+16n(3 + 2n)A0 - 16A 2 ] + 8 n(n + 2)A} / (32P'(0) 3 } (32) 

This form of the Ricci scalar shows that the singular points r smg of a given solution (15)-(16) 
of the EGB equations either satisfy lim r _ >rsjns 0(r) = ±oo or lim r ^ rs , ng 0(r) = 0 O , 0 o being 
the stationary point of P(ip). If /x = 0 then ip(r) = constant and the horizon is trivially 
found. In the k = 0, /a ^ 0 case there will be a horizon only if -0 = 0, which requires that /j 
and A have opposite signs. The horizon will be at 

r h = (~(n + l)/i/A) 1/(n+1) (K = 0,g^0) (33) 

If /i ^ 0, k = ±1, there is a horizon at every point where 

sgnf'ip) — k and P(i/)) = — (k — ±1 ,n ^ 0 ), (34) 

a 

For later convenience, we re-write (16) as 

(35) 

where 


a, =- L -7-i±./i + !W y- 1 )(’ i 

a(n — l)(n — 2) l y n(n + 1) 

Note that, for /i/cx > 0 the condition f — k — r 2 ^ > 0 reduces to 



-0 < 0 or 0 < 0, —10| 2 < P(0) (if k — 1) 
a 

0 < 0 (if k — 0 ) 

0 < 0 and — 10| “ 2 “ > P(0) (if k — —1) 
a 

whereas for /i/a < 0 , f — k — r 2 0 > 0 is equivalent to 


0 < 0 or 0 < 0, P(0) < —10| 2 (if k — 1) 

a 

0 < 0 (if k — 0 ) 

0 < 0 and P(0) > —10| “ S - (if k — —1) 


(36) 


(37) 

(38) 

(39) 


(40) 

(41) 

(42) 



7 


We label solutions with a three digit number in the form a.b.c, with a, b and c labeling the 
distinct ranges of values for a, A and p respectively. Only the 1.1.c cases (positive a and A) 
will be analyzed in full detail, since the other cases are trivial variations of this one. A plot 
of P and fjb\ip\( n+1 ')/ 2 /a vs ip is given in each case; the positive (negative) ip intersections of 
these curves are k — 1 (k — —1) horizons, ip = 0 being the horizon when k — 0. If p/a > 0 
(< 0), the portion of P above (below) the ip axis gives the two solution branches ipipr) of 
(16). Note from (16) that r can extends to infinity only if P has real roots, and that there 
is a singularity at r sing if either lim r _ vrsjn9 ippr) = ±oo or lim r ^ rsing ip(r) = ip Q , P’(ip 0 ) = 0 . 
Equations (37)-(39) ((40)-(42)) are used to find the / > 0 region of interest when p/a > 0 

(< o). 

Case 1 : a > 0 
Case 1.1: A > 0 

In this case P has two real roots Ai < 0 < A 2 , |A 2 | < |Ai|. If p > 0 (35) has two solutions 
ipipr), i = 1,2, with r extending to infinity, and linv^^ ipipr) = A t (Figure l.a). For ipipr), 
as r goes from infinity down to zero, ip runs from Ax to —oo, where, according to our 
previous analysis, there is a curvature singularity. Similarly, for ip 2 pr), as r goes down to 
zero, ip runs from A 2 to +oo, where there is a curvature singularity. Three qualitatively 
different p > 0 cases are plotted in Figure l.b. The positive ip intersections with P give 
horizons in the n = 1 case (eq.(34)), the negative ip intersections with P give horizons in 
the k — —1 case, and ip = 0 is the n = 0 horizon. Note from (34) that some of the drawn 
k = ±1 horizons (curve intersections) may be missing in the special case n = 3. Note also 
that, as A —> 0, A 2 —> 0 and some horizons move to infinity. 

Case 1.1.i: large positive p (figure l.a and curve (i) of figure l.b) 

In view of eqns (37)-(39) the ip 2 branch never gives / > 0, whereas, for any n, the ip\ branch 
gives a spacetime with r sing — 0 < r < oo (naked singularity). 

Case 1.1.ii: intermediate positive p (figure l.a and curve (ii) of figure l.b) 

The analysis for the ipi branch is as in case 1.1.i. The ip 2 branch gives / > 0 for k — 1, 
case in which the spacetime has two horizons and no singularities, rp on < r < rh or2 ■ As 
A —> 0 + , Tf lor2 —> oo. If n — 3 one of the intersections of P with curve (ii) may be absent, 
and r sing = 0 < r < r hor . 

Case 1.1.iii: small positive p (figure l.a and curve (iii) of figure l.b) 

The analysis for the ip 2 branch is as in case 1.1.ii. For ipi and k — 0,1, r svng = 0 < r < oo 
(naked singularity), whereas for n = —1 there are two / > 0 regions, one for which 
'r S ing — 0 < r < rh or i (naked singularity), the other satisfying rp or2 < r < oo . The first 
region may be missing if n — 3. 

Case 1.1.iv: p = 0 

We obtain cosmological, non-singular solutions /(r) = 1 — r 2 A 2 , with 0 < r < A/ 1 ^ 2 , 
f(r ) — k — r 2 Ai, k, — 0,1 and 0 < r, and /(r) = —1 — r 2 Ai, (k = —1), r > IAxI -1 / 2 . 

If p < 0 (35) has two solutions ipipr), i = 1,2, with r extending to infinity, and 
linp^oo ipipr) = A*. There is a minimum value r = r sing defined by ip\ (r smg ) = ip 2 pr smg ) = ip 0 
pipe = (Ai + A 2 )/ 2 ), this point is singular in view of eq. (32) because P\ip 0 ) = 0 . As r grows 
from r sing to infinity, ipi goes from ip Q to A* (figure 2 .a). 

Case 1.1.v: small negative p 

The ipi branch has, for n = —1, a horizon that hides a singularity, / > 0 if pr sing <)rh or < 
r < oo, whereas for n — 0,1 there is a naked singularity, r sing < r < oo. The ip 2 branch gives 




FIG. 1: Cases 1.1.i to 1 . 1 .iii: a) the two branches = 1,2 of equation (35) in the case 

li/a >0. V’i —► Aj as r —> oo, V’l (V^) tends to —oo (+oo) as r —■> 0 + . b) Plots of fj,\^ n+1 ^ 2 /a for 
(i) large (ii) intermediate and (iii) small positive n/a. 


no / > 0 solution for k = — 1 , whereas for k — 0,1 gives a spacetime with r sing < r < Thor- 

Case 1.1.vi: large negative // 

For any k, the ipi branch gives a space-time with a naked singularity, / > 0 for r smg < r < oo. 
For any k, there is a horizon in the ip 2 branch, and / > 0 for r sing < r < Thor- 



FIG. 2: Cases 1.1.v and 1.1.vi: a) the two branches = 1,2 of equation (35) in the case 

n/a. < 0, V’i goes from ip 0 to Aj as r goes from r S i ng to oo. b) Plots of P and i 2 \ip\^ n+1 ^ 2 /a vs. i/i 
for (v) small and (vi) large negative fi/a. 

Case 1.2: — n(n +1)/ (4a(n — l)(n — 2)) < A < 0 

In this case P has two real roots A x < A 2 < 0. Six cases of /i values should be distinguished: 
1 . 2 .i large positive, 1 . 2 .ii small positive, 1 . 2 .iii null, 1 . 2 .iv small negative, 1 . 2 .v intermediate 
negative and 1.2.vi large negative. These different cases are reepresented in figure 3 below. 
Case 1.3: A < — n(n + 1)/ (4 a{n — 1 )(n — 2)) 

P has complex roots, n/a must be positive, and there is a maximum value of r (corresponding 
to 'ip — i \) 0 ) which is singular. Three ranges of /x values should be distinguished: 1.3.i 
large positive, 1.3.ii intermediate positive and 1.3.iii small positive. These are illustrated in 
figure 4. 
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FIG. 3: Cases 1.2.i to 1.2.vi: a) the two branches 0j(r),z = 1,2 of equation (35) in the case 
n/a > 0 together with the fj,\ip\( n+1 ^ 2 / a curve for (i) large and (ii) small positive values of fi/a. 
b) ipi(r),i = 1,2 in the case /.i/a < 0 together with the n\^ n+1 ^ 2 /a curve for (iv) small (v) 
intermediate and (vi) large negative values of n/a. 



FIG. 4: Cases 1.3.i to 1.3.iii: P and fi\ip\( n+1 ^ 2 / a for (i) large, (ii) intermediate and (iii) small 
positive values of [ija. P has no real roots, as r grows from r S i ng , ipi(r) (02(r)) moves to the left 
(right) of 0 O . 


V. STABILITY OF MAXIMALLY SYMMETRIC STATIC SOLUTIONS 


The stability of the solutions (1)-(13)-(16) of the EGB vacuum equation can be analyzed 
using the “S-deformation” approach [16]: consider the operator 


A : = 


d 2 

dr* 2 




(43) 


acting on smooth functions defined on / = {r*|rj[‘ < r* < r * 2 }, the regular, / > 0 region 
(note that it is possible that r* = ±cx)). E in (28) is greater than or equal to the lower 
bound of (0, A0)/(0, 0), 0 smooth of compact support on /. However, for any such 0, given 
a smooth S, 



(|D0| 2 + U|0| 2 ) dr*, 


(0, A<j>) 


(44) 
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where 


and the “deformed potential” V is 


D = ^ + s ' 


V = V+ - S 2 

dr 


(45) 


(46) 


If an S function is found such that V > 0 on /, the stability of the solution is guaranteed, 
as follows from (44). Note from (31) that the choice 


S = -A ln(A) 
dr 


gives V = q, then 


Defining 


{<j>,A<j>)= / \D(j)\dr* + 


1 7*1 


/ 


dr. 


H = 


9 

r q 


/(2k- 7 ) 

the expectation value of A can be conveniently written as 

(<p,A(j))= [ \D(p\ 2 dr* + ( 2 k- 7 ) [ 


2 H 


■ dr, 


(47) 


(48) 


(49) 


(50) 


Note that neither H nor D depend on 7 . This factorization of the “deformed potential” q 
is the one referred to in Section I, and is crucial to arrive at the stability criterion below. 

If the Riemannian manifold E” is compact without boundary, applying Stokes’s theorem to 


0 < / (VW fc - n\/ j h ik ) ( Vih jk - kV jh ik ) 


(51) 


'S2 


and using the TT condition of h tJ together with Rijki = n(gik9ji — 9jk9u) we arrive at 

7 < — n 2 n =>• 2k, — 7 > K 2 n + 2k > 0 (52) 

for n > 3. Then from (50) we conclude that H > 0 on / implies stability. Now suppose 
H < 0 at some point in /, then a test (j) can be found such that 


2 H 


dr < 0 . 


(53) 


' 7*1 


The ‘kinetic” piece of (50) may certainly be larger than the absolute value of the integral in 
(53), but (50) will be negative for sufficiently high harmonics. We conclude that a solution 
is stable if and only if H > 0 on I. 

Using (35) and its first r derivatives, and introducing = (Ai + A 2 )/2 and A = (A 2 — 
Ai)/2, a simple expression for H in terms of if is found which is even in x = (ip — if 0 )/ A: 


H = 


(n - 3) (n - 5) x 4 + 2 (n + 1) (2 n - 3) x 2 - (n + l ) 2 
2 (n — 2) x 2 ( x 2 (n — 3) + (n + 1)) 


(54) 
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Also, if fx 7 ^ 0 , 



n +1 


a 




1 

n+1 


n(n — l)(n — 2 )A 2 ^ "+ 1 

4 J 


x la ; 2 — 1 


l/(n+l) 


H 


x* 


1 ) 


dx (55) 


An immediate consequence of the stability criterion above and (54) is that the EGB 
cosmologies are all stable against tensor perturbations, since, for (i — 0, ^( r ) = A, : , then 
x = ±1 and H — 1. 

Note that the cases n = 3,4, 5 of (54) are special: 


H(n= 3 ) 
H(n= 4 ) 
H(n= 5 ) 


3 X 2 — 2 

X 2 

-a : 4 + 50a: 2 - 25 
4a : 2 (a : 2 + 5) 
7a: 2 -3 
a : 2 (a : 2 + 3) 


(56) 

(57) 

(58) 


A. Stability analysis 


When P has real roots, a: is real in (54)-(58). Furthermore 


P(n=3) >0 iff | x | > 


~ 0.82 


H (n=4) > 0 iff 0.71 ~ a/15 - VTO < |a:| < \/l5 + VTo ~ 7.03 


H(n= 5 ) > 0 iff |a:| > 


~ 0.65 


H, 


(n> 5 ) 


> 0 iff |a:| > 


(n — l)(2n — 3) 
(n — 3)(n — 5) 


(n — 3) (n — 5) (n + 1)" 


1 + 


(n - l) 2 ( 2 n - 3 ) 2 


The r.h.s. of (62) decreases to a/—: 2 + — 0.49 as n grows from n = 5. 


(59) 

(60) 
(61) 



Cases 1.1.i to 1.1.iii: 

All these solutions are stable if n ^ 4. The stability follows from (59), (60) and (62) above, 
which show that H > 0 if |a:| > 1 (i.e., if} > A 2 or -0 < A x ) . The case n — 4 is special, as 
follows from (60) and was already noticed for k — 1 and A = 0 in [15]. We now analyze the 
stability of every cosmological and black hole n — 4 solution found in cases 1.1.i through 
1 . 1 .iii: 

Cases l.l.ii-iii, branch, n — 1: this black hole solution has r S i ng < rh ori < r < rhor 2 and 
will be stable as long as Xh ori < \/l5 + VTO, i.e., for large enough /i. Note that this is also 
true for A = 0 (A 2 = 0), the low // instability being the one found in [15]. 

Case l.i.iii, ipi branch, k — — 1: there are two / > 0 regions, and one of them gives a black 
hole solution, which will be stable as long as Xh or < \/l5 + \/TC)- Contrast this condition to 
that obtained before, k — —1 black holes require small (i to be stable. 
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Case 1.1.iv 

As explained above, H — 1 in this case. These cosmological solutions are all stable. 

Cases 1.1.v to 1.1.vi 

From (59)-(62) the n = —1, if )i branch black hole in case 1.1.v will be stable for \fi/a\ small 
enough. 

The analysis of the remaining cases can be readily done as in the previous cases and is 
left to the reader. 


VI. CONCLUSIONS 

We proposed a classification scheme for static solutions to the Einstein Gauss Bonnet 
gravity of the form E” x R + , E™ an n— manifold of constant curvature k, and studied 
their linear stability under tensor mode perturbations. We found an explicit form of the 
potential of the Schrodinger like equation governing the time evolution of the perturbation, 
and studied its spectrum using the A—deformation approach. An S'—deformed Schrodinger 
potential was found that conveniently factors out the eigenvalue of the laplacian on E™ 
associated with the perturbation, allowing a definite classification of every space-time into 
stable or unstable. Preliminary results indicate that this feature of tensor perturbation is 
shared by vector perturbations. The scalar case is still under investigation [17]. Cosmological 
solutions and a variety of Euclidean or dS black holes with positive curvature horizons are 
shown to be stable in space times of dimension d ^ 6 with a positive values of a -the Gauss 
Bonnet term coupling-. In six dimensions, these black holes are stable only if their masses 
are above a critical value. Black holes with negative curvature horizons are found in any 
dimensions which are stable only if their masses are below a critical value (see [22] for a 
thermodynamic instability of black holes in EGB gravity). 

The stability of these space times under vector and scalar perturbations is currently being 
studied. 
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APPENDIX A: LINEARIZATION FORMULAS 


The first order variation of the Riemann tensor, Ricci tensor and Ricci scalar under (4) 
can be obtained after a long calculation using (4)-(11). These are: 


6R fi tj = 


f 


2 / 4 fg 


hj 


SR ti rj = ~ + 


r 


2 g \ 4 /.9 2 gr 

S R ti jk = (v fc v - V j h k ) 

r 


</> W 


SR ri tj = 
5R r P = 


2 fr 4 p 


4 g 1 2 rg 


+ 


2 / 

Jf 


5Rri jk = (V V - W j h k ) 


2 g 


2 r 2 


hi j 


SRi/ k = JJ (Vihj k - vfhf) 


= — 

2 ^ 




(V,/;/- - V,/>/') 

|§)+ Ij] “ S ‘R + W ~ 4U) 

+ 2^ + VW - V,-V V) 


(Al) 

(A 2 ) 

(A3) 

(A4) 

(AS) 
(A 6 ) 

(AT) 
(AS) 


(A9) 


the other components of 5R a b cd being zero. The nonzero components of the Ricci tensor 
then are 


5Ri j = 

4+u 

f sf 

n f \ 

0 " 

+ -^7 (2 K — 7 ) 


2 / V 

V 4 g 2 

2rg 4 fg) 

"2 g 

2 r 2 V 


h j 

IL l 1 


Finally, 

SR = 0 

From these equations and (14) follows (20). 


(A10) 


(All) 
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